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Abstract With the assumption of a linear-dashpot interac- 
tion force, we compute the coefficient of restitution as a 
function of the elastic and dissipative material constants, k 
and 7 by integrating Newton's equation of motion for an iso- 
lated pair of colliding particles. It turns out that an expres- 
sion e(k, y) widely used in the literature is even qualitatively 
incorrect due to an incorrect definition of the duration of the 
collision. 
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1 Introduction 

The (normal) coefficient of restitution of colliding spherical 
particles relates the post-collisional velocities and v'- and 
the corresponding pre-collisional velocities v, and v,-, 



/ l+£ w \ 1 

V,- =V,--^— [(Vf-V/J-eyJey 



v ;= v ; 



where 



2 

1 + e 



(1) 



[(▼i-Vj)-eyK 



r, r, 



(2) 



This research was supported by German Science Foundation (Grant 
PO 472/15-1). 

Thomas Schwager 
Charite 

Augustenburger Platz 
10439 Berlin 
Germany 

E-mail:thomas.schwager@charite.de 

Thorsten Poschel 
Charite 

Augustenburger Platz 
10439 Berlin 
Germany 

E-mail: thorsten.poeschel@charite.de 



is the unit vector at the instance of the collision. For simplic- 
ity, we consider particles of identical mass - the generaliza- 
tion to unequal masses is straightforward. 

The coefficient of restitution is the most fundamental 
quantity in the theory of dilute granular systems. This co- 
efficient together with the assumption of molecular chaos 
is the foundation of the Kinetic Theory of granular gases, 
based on the Boltzmann or Enskog kinetic equation, see, 
e-g-! EL 0; HI] an d many references therein. On the other 
hand, the coefficient of restitution is also essential for event- 
driven Molecular Dynamics simulations which allow for a 
significant speed-up as compared to traditional force-based 
simulations, see, e.g., 19D for a detailed discussion. The most 
important precondition of Eq. (QJ is the assumption of ex- 
clusively binary collision, that is, the system is assumed to 
be dilute enough such that multiple-particle contacts can be 
neglected. The latter condition seems to be rather restrictive, 
however, it was shown in many simulations that event-driven 
simulations are applicable up to rather high density ]6[]. 

The coefficient of restitution is, however, not a funda- 
mental material or particle property. Instead, the particle in- 
teraction forces and Newton's equation of motion govern the 
dynamics of a mechanical many-particle system. Therefore, 
if we wish to use the concept of the coefficient of restitu- 
tion and the simple equation (QJ to compute the dynamics 
of a granular system, we have to assure that Eq. ([]]) yields 
the same post-collisional velocities as Newton's equation of 
motion would do. 



2 Coefficient of restitution and Newton's equation of 
motion for frictionless spherical particles 

Let us consider the relation between the coefficient of resti- 
tution and the solution of Newton's equation of motion. As- 
sume two approaching particles at velocities vj and V2 come 
into contact at time f = 0. From this instant on they deform 
one another until the contact is lost at time t =t c . We shall 
consider the end of the collision at time t c separately below. 
For sufficiently short collisions we can treat the unity vector 
e as a constant. Then we can describe the collision by the 
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mutual deformation 



with initial conditions 
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where ri (f ) and r 2 (?) are the time-dependent positions of the 
particles and R is their radius. The interaction force between 
the particles is model specific, e.g. I^ fllll . and, in general, a 
function of the material parameters, the particle masses and 
radii, the deformation £, and the deformation rate | . Its time 
dependence is expressed via F = F {t (t), | (t)\ Newton's 
equation of motion for the colliding particles reads 



(4) 



where v is the normal component of the initial relative ve- 
locity, 



[vi(0)-V2(0)]-[ri(0)-r 2 (0)] 
2R 
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The coefficient of restitution follows from the solution of Eq. 



Ik), 
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where t c is the duration of the collision. 

Thus, the coefficient of restitution due to a specified in- 
teraction force law can be determined by solving the equa- 
tion of motion |[l0l;[TTI;[l3ll. In general, £ is a function of the 
material parameters, the particle masses and radii and the 
impact velocity v. 



3 Linear dashpot model 

Consider the linear dashpot force, 



(7) 



This force is problematic as a particle interaction model since 
particles made of a linear-elastic material do not reveal a lin- 
ear repulsive force, neither in tree dimensions [4] nor in 2D 
JH]. The subsequent analysis can be also performed for more 
realistic force models, such as viscoelastic forces in 3d |[2j] 
and 2d lfl2ll which will be published elsewhere. 

The linear dashpot model is of interest here because this 
force leads to a constant coefficient of restitution, e.g. fTHl . 
that is, e does not depend on the impact velocity but only 
on the material parameters k and y. The constant coefficient 
of restitution in turn is the preferred model in both the Ki- 
netic Theory of granular gases and also in event-driven sim- 
ulations of granular matter. In contrast, more realistic force 
models lead to an impact-velocity dependent coefficient of 
restitution @; M, d S QJ] . 

The equation of motion corresponding to Eq. (0 
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has two solutions. With the abbreviations 

co^-; = 1; «^A 2 -j3 2 . 
m m v 
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we can write the solution for the case of low damping (/3 < 
(Oq) as 



§(0 = ---* 



CO 



e sin COt , 



while for the case of high damping (/3 > COq) we have 

£(f) = sinflf, 

with 
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To determine the coefficient of restitution according to Eq. 
© we need the duration of the collision t c . A natural choice 
seems to be 



«(£)=0; 



tf>0, 
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that is, the collision is finished when | r, — Tj \ = 2R. With this 
condition we obtain t® = n/co and the coefficient of restitu- 
tion 
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for the case of low damping which is widely used, e.g. in 
flDl and many others. 

Unfortunately, condition dl4l) cannot be correct since at 
time t^l = n/co the interaction force is negative, 



\coJ CO 
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see Fig. Q] which contradicts the assumption of exclusively 
repulsive interaction between granular particles. Moreover, 
for the case of high damping Eq. (114b has no real solution, 
that is, the particles stick together after a head-on collision 
(see Fig.|2]l. One can call this situation dissipative capture. 

To avoid these artifacts, instead of Eq. (114b . we now ex- 
plicitely take into account the condition of no attraction. 
That means, instead of Eq. © we have to use the force 
F* = max(0, F), with F given by Eq. (0. This expression 
as is in common use in (force-based) Molecular Dynamics 
simulations, assures that only repulsive forces act. It is clear 
that once the force becomes zero, the collision is finished, 



|(fc) = 0; t c >0; |(f c )<0, 



(17) 



that is, the collision is complete when the repulsive force be- 
tween the particles vanishes and the surfaces of the particles 
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Fig. 1 Sketch of a particle collision. The upper part of the figure visualizes the condition <14> for the end of the collision. As an artifact of the 
model, there appear attracting forces towards the end of the collision (boxed notation). Lower part: The condition U7\ for the end of the collision 
avoids this artifact. 



separate from one another. This condition takes, thus, into 
account that Eq. © applies to particles in contact. 

For the case of low damping Eqs. ( 111b and ( 117| l yield 
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One has to take care to select the correct branch of the arc 
tangent for j8 < C0q/ \fl (or j8 < ft)) and for /3 > (Oq/s/1. The 
solution reads 
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For the coefficient of restitution we find 
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regardless of the branch of the solution of Eq. ( 1181) . Together 
with the solution for the case of high damping which can 
be obtained by straight-forward computations we obtain the 
coefficient of restitution: 
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In accordance with physical intuition we have > £® since 
the result given in Eq. (1211) does not suffer from the artifact 
of unphysical attractive forces implied in Eq. dl5l) . The un- 
physical dissipative capture does not occur - even for large 
damping constant the particles separate eventually. For very 
high damping the coefficient of restitution can be approxi- 
mated as 



ft). 



(22) 



£ d~^2 for /3>ftJo 

Fig. shows both the (wrong) coefficient of restitution 
as derived from criterion ( 114b (dashed line) and the coef- 
ficient of restitution derived from the correct criterion ( 117b 
(full line). As the solution Eq. (12TT) only depends on the ratio 
a = /3 /(Ho we draw the coefficient of restitution as a func- 
tion of this parameter. Figure [2] shows that the unphysical 
attraction indeed yields a too low coefficient of restitution. 
The inset shows a logarithmic plot of the same curve show- 
ing that there is no capture if the correct criterion is applied. 



4 Conclusion 

Assuming a linear dashpot interaction force we computed £ 
as a function of the elastic and dissipative material proper- 
ties by integrating Newton's equation of motion for a pair of 
colliding particles. The condition for the duration t c of the 
collision <Sj (r!) = 0, that is, the interaction stops when the 
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Fig. 2 Coefficient of restitution as function of a = P/coq. The dashed 
curve is the result for the wrong termination condition 0. Obvi- 

ously, for high damping a > 1 we have dissipative capture. The inlay 
shows ln£ as a function of a for the correct termination condition. 
There is no dissipative capture even for very high damping. 

particles have the distance |r,- — r ; -| = 2R, is widely used in 
the literature, however, it leads to erroneous attracting forces 
between the particles and, thus, to wrong values for the co- 
efficient of restitution. In extreme cases, the result may be- 
come even qualitatively incorrect since despite the interac- 
tion force being purely repulsive, the particles may agglom- 
erate. In this state one would need to supply energy to sepa- 
rate the particles from one another. 

To account for the fact that the interaction force can never 
become attractive, we use the condition fc, (t c ) = which 
avoids the mentioned artifacts. This condition takes into ac- 
count that the collision may be completed even before ^ = 0, 
that is, the surfaces of the particles lose contact slightly be- 
fore the distance of their centers exceeds the sum of their 
radii. Thus, the deformation of the particles may last longer 
than the time of contact and the particles gradually recover 
their spherical shape only after they lost contact. 

The resulting coefficient of restitution is always larger 
than the value based on % = 0. In agreement with the 
repulsive character of the interaction force, the coefficient of 
restitution is always well defined and non-zero. 
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